Stochastic resonance in soft matter systems: combined effects of static and dynamic disorder 



Matjaz PercQ Marko Gosak[j] and Samo Kraljl 

Department of Physics, Faculty of Natural Sciences and Mathematics, University of 
Maribor, Koroska cesta 160, SI-2000 Maribor, Slovenia 
* Jozef Stefan Institute, P.O. Box 3000, SI-1000 Ljubljana, Slovenia 

We study the impact of static and dynamic disorder on the phenomenon of stochastic resonance (SR) in a 
representative soft matter system. Due to their extreme susceptibility to weak perturbations soft matter systems 
appear to be excellent candidates for the observation of SR. Indeed, we derive generic SR equations from a 
polymer stabilized ferroelectric liquid crystal (LC) cell, which is a typical soft matter representative constituting 
one of the basic components in several electro-optic applications. We generalize these equations further in order 
to study an even broader class of qualitatively different systems, especially disclosing the influence of differ- 
ent types of static disorder and interaction ranges amongst LC molecules on the SR response. We determine 
the required conditions for the observation of SR in the examined system, and moreover, reveal that a random 
field type static disorder yields qualitatively different responses with respect to random dilution, random bond 
and spin glass universality classes. In particular, while the latter three decrease the level of dynamic disorder 
(Gaussian noise) warranting the optimal response, the former evokes exactly the opposite effect, hence increas- 
ing the optimal noise level that is needed to resonantly fine-tune the system's response in accordance with the 
weak deterministic electric field. These observations are shown to be independent of the system size and range 
of interactions, thus implying their general validity and potentially wide applicability also within other similar 
settings. We argue that soft matter systems might be particularly adequate as a base for different SR-based 
sensitive detectors and thus potent candidates for additional theoretical as well as experimental research in the 
presently outlined direction. 
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I. INTRODUCTION 

When noise is introduced to nonlinear systems one 
can observe a variety of interesting and counterintuitive 
phenomena. 1,2,3 Perhaps the most prominent is the phe- 
nomenon of SR where an appropriate intensity of noise evokes 
the best correlation between a weak deterministic stimulus 
and the response of a nonlinear system, 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 1 ¥ 2 - 13 This 
contradicts intuitive reasoning that suggests noise can only 
act destructive on system performance. To observe the phe- 
nomenon of SR a nonlinear systems must contain the follow- 
ing three basic ingredients: i) an energetic activation barrier 
{i.e. a threshold), ii) a weak coherent input (i.e. periodic ex- 
ternal forcing), and iii) a source of noise. These conditions 
are very general indeed, enabling the observation of SR in the 
most diverse circumstances ranging from neuronal and brain 
functioning 14 ' 15 ! 16,17 ' 18 ! 19 ! 20 to the recurring occurrences of ice 
ages, as comprehensively reviewed in the past i 21 ' 22,23 ' 24 The 
SR was first verified experimentally on an electronic circuit 
with a binary switch, 25 but latter on many additional exper- 
imental findings were reported, as for example in a bistable 
non-linear optical device. 26 

Directly linked to the subject of the present work are stud- 
ies focusing on noisy bistable oscillators* 2 - 7 ^ as well as stud- 
ies that consider coupled bistable systems 2 9 ' 30 ' 31 ' 32 ' 33 within 
different settings of SR. The latter setup served also as a basis 
for the discovery of an interesting effect of noise on an ensem- 
ble of coupled bistable oscillators; namely the so-called array 
enhanced SRi^ 5 . Thereby, Lindner et al. demonstrated that 
a linear coupling, combined with noise and a weak periodic 
signal, can enhance synchronization and global organization 
in a chain of overdamped nonlinear oscillators. Notably, in- 



vestigations of bistable noisy systems remain vibrant to date, 
and in the last few years several theoretical 38 ' 39 ' 40 as well as 
experimenta l 36,37 studies have been performed that report on 
interesting new aspects of SR in such systems. Furthermore, 
noise-supported signal prop ag ation 4 1 ' 42 ' 43 and the SR of col- 
lective variables with very large gains 46 in arrays of bistable 
systems have been studied as well. Recently, the scope of SR 
and related phenomena in coupled oscillators shifted also to 
ensembles characterized with complex interaction topologies, 
as constituted by small-world or scale-free networks. 44 ' 45,47 

Since soft matter systems 48 evolve both in time and space, 
and are in addition to that extremely susceptible to external 
perturbations, it is rather remarkable that they have avoided 
being the subject of investigations related to SR in the past. 
In particular so since during the last two decades the phe- 
nomenon of SR has been reported in the most diverse circum- 
stances, ranging from semiconductor devices, 49 magnetoelas- 
tic ribbons, 50 superconducting resonators, 51 intercellular cal- 
cium wave dynamics, 52 SQUID based detectors for weak 
magnetic fields, 54 to sociological models. 53 Moreover, the SR 
effect has also been reported in chemical reactions, where it 
has been shown that optimal noise intensities are able to en- 
hance weak periodi c 56 ' 57 - 58 - 59 as well as aperiodic 60 signals, 
or even support traveling waves in a spatially extended chem- 
ical medium. 61 Recently, the popular nanoscale systems be- 
came associated with the SR as well. In particular, it has been 
shown that nanomechanical oscillators in a dynamic bistable 
state exhibit a more controllable switching in the presence of 
noised 5 - 

At present, we aim to extend the scope of SR also to soft 
matter systems, considering a polymer stabilized ferroelec- 
tric (PSF) LC cell as a typical soft matter representative. Al- 
though the phenomenon of SR is most commonly associated 
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with dynamic disorder or noise, soft matter systems, and in- 
deed real-life systems in general, are in addition characterized 
by a certain degree of static disorder. This fact has recently 
been addressed by Tessone et al.r^ who succeeded showing 
that static or quenched disorder can give rise to SR as well. 
Here we examine the impacts of both static and dynamic dis- 
order on the phenomenon of SR, and demonstrate under which 
conditions it can be observed in a PSFLC cell. More precisely, 
we consider different origins of static disorder and show that a 
random field type static disorder yields qualitatively different 
responses with respect to random dilution, random bond and 
spin glass universality classes. In particular, while there al- 
ways exists an optimal level of dynamic disorder, warranting 
the best correlation between a weak periodic electric field and 
the response of the PSFLC cell, the considered types of static 
disorder may additionally enhance or deteriorate the SR. Ran- 
dom dilution, random bond and spin glass universality classes 
decrease the level of dynamic disorder warranting the opti- 
mal response, whereas the random field evokes exactly the 
opposite effect, hence increasing the optimal level of dynamic 
disorder that is needed to resonantly fine-tune the system's 
response in accordance with the weak deterministic electric 
field. We additionally test these findings on their robustness 
with respect to the system size and range of interactions, and 
reveal that they are largely independent of such particularities. 
We discuss that, due to their extreme susceptibility to exter- 
nal perturbations and wide applicability, soft matter systems 
might represent an interesting environment for the develop- 
ment of SR-based sensitive detectors, and moreover, that ad- 
ditional theoretical as well as experimental research efforts in 
the direction outlined presently seem justified to realize these 
potentials. 

The remainder of this paper is organized as follows. In Sec. 

II we present a generic dynamic equation giving rise to the SR 
phenomenon and the simulation procedure to solve it. Section 

III is devoted to the derivation of governing equations of the 
PSFLC cell, whereas in Sec. IV we present the results. In par- 
ticular, we estimate analytically the influence of static disorder 
on the system's configuration, and numerically investigate the 
combined influence of static and dynamic disorder on the SR. 
The final section features the summary of presented results 
and outlines possible applications of our findings. 

II. STOCHASTIC RESONANCE: GENERIC EQUATIONS 

We consider a system of N coupled bistable overdamped 
oscillators, governed by the Langevin equations of the form 

2* = Si - si + J^2( Sj - Sl )+E + V2D&(t), (1) 
j 

where Si describes the state of the i-th oscillator, t is the di- 
mensionless time, the sum index j runs over the oscillators sj 
that are coupled with Si, J measures the coupling strength 
amongst the oscillators, E = Eo cos(wt) is a weak peri- 
odic field oscillating with the frequency u>, whereas 2D is the 
variance of Gaussian noise with zero mean and autocorrela- 
tion (£,i(t)^j(k)) = SijStk- Equation ((TJ arguably provides a 



paradigmatic setup for different scenarios of SR, either via the 
classic setup, or via variations of system size^ or diversity. 62 
In the absence of deterministic forcing and noise each oscilla- 
tor is characterized by a double-well potential with minima at 
Si = so = ±1, whereas a periodic field E with finite ampli- 
tude Eq > modulates the double-well potential, potentially 
breaking the reflection symmetry of the system. In the regime 
of our interest Eq is too weak to let the oscillators roll period- 
ically from one potential well into the other. Instead, the min- 
ima are alternatively and asymmetrically pushed up and down, 
periodically raising and lowering the potential barrier. By in- 
creasing D the noise-induced hopping between the potential 
wells can become synchronized with the frequency u of the 
weak external field, thus exhibiting the SR phenomenon. 

To evaluate responses of the noisy system to the periodic 
field E in dependence on D, we calculate the Fourier coeffi- 
cients for the mean field 

i 

Noteworthy, the Fourier coefficients are proportional to the 
(square of the) spectral power amplification, 64 which is of- 
ten used as a measure for the quantification of SR. Here, the 
Fourier coefficients Q are calculated according to the equa- 
tions 

,27m/u 

Q si „ = — / S(t)sm(cjt)dt, 
r2*n/u> 

Qcos = — / S(t) cos(wt)di, (3) 

™ Jo 

Q = \foL + oL, 

where to is the frequency of the weak periodic field E and n 
is the number of oscillation periods used. 



III. STOCHASTIC RESONANCE IN THE PSFLC CELL 

In order to demonstrate an appropriate soft matter system 
that can potentially exhibit SR, we start by considering a ther- 
motropic bistable surface stabilized ferroelectric (SSF) LC 
cellp£ The ferroelectric LC, exhibiting both translational and 
orientational (quasi) long-range order, is confined in a plan- 
parallel cell as shown in Fig. Q] The translational order is 
characterized by smectic layers that are stacked along the z- 
axis, and the normals of the cell walls are parallel with the 
x-axis of the system. The orientation ordering is described 
by the director field n, yielding an average local orientation 
of LC molecules, where |n| = 1. The tilt (also referred to 
as the cone) angle 8 = arccos(n • v) between the smectic 
layer normals n and v distinguishes between the smectic A 
(SmA) and chiral smectic C (SmC*) ordering. The SmA or- 
dering, in which (9) = 0, is commonly realized above the 
critical temperature T c , whereby (■) stands for the spatial aver- 
age. Below T c the LC enters into a helicoidal SmC* ordering. 
In a bulk sample (0) > the direction of the tilt precesses 



3 



j£jr'Jr_b i JrirJr 

_b' J? _2 _b J: Jr jJr _b 
1- 1- JLk 2 1- 1- 1- 
t* b r b:bir_bJr 
_b Jr _b _b ir Jr ± _b 
i: i: _b _b Jb _b _b _b 



FIG. 1: The SSFLC cell and the geometry of the problem. LC 
molecules are collected in smectic layers that are stacked along the 
z-axis. For d < p the molecules are either tilted along the positive 
or the negative y-axis. Consequently P is aligned along the negative 
or the positive i-axis, respectively. The nail representation is used 
in order to show that LC molecules are tilted in the (y, z) plane. Ar- 
rows attached to the 'nails' mark the local polarization of molecules, 
whereas muddy gray lines in the background roughly indicate the 
structure of the polymer network. The weak periodic electric field E 
is introduced along the x-axis. 



as one goes from one smectic layer to another. The result- 
ing helicoidal structure is characterized by the wave number 



q = 2n/p, where p stands for the pitch of the helix. Tilt- 
ing of the molecules gives rise to the in-plane spontaneous 
polarization P, which is in general perpendicular to n and v, 
whereby due to symmetry reasons it holds P = Pv x n. 6b 
Therefore, in the bulk SmC* phase P also forms a helix and 
the macroscopic polarization of the system equals zero. The 
cell confinement can suppress the helicoidal structure if the 
thickness d of the cell is less or comparable to p. The re- 
sulting unwound SmC structure is on average homogeneously 
tilted along a single direction. In the case of isotropic planar 
anchoring (i.e. molecules tend to lie in the confining plane 
within which all directions are equivalent), the molecules tilt 
either along the positive or the negative y-axis. We henceforth 
refer to these configurations as Up and Down states, which 
are separated by an energy barrier. Importantly, such a system 
thus exhibits bistability. The resulting polarization is aligned 
either along the positive or the negative x-axis, respectively. 
If a periodic external electric field is applied along the a>axis, 
it alternately favors the Up and Down configuration. If also 
a noisy component is present, then all the three essential con- 
ditions for the SR phenomenon listed in the Introduction are 
fulfilled. 

We proceed by describing the bistable equilibrium struc- 
ture in the unwound SmC phase quantitatively in terms of 
the tilt angle 9 and the polarization P. We consider the sys- 
tem depicted in Fig. Q] where p < d. The molecules are ei- 
ther in the Up or Down state, i.e. n = (0, | sin0|,cos0) and 
P = (=f|P|, 0, 0). We further apply an external electric field 
E = Eo(cos(uit), 0, 0). We expand the corresponding free 
energy density / in terms of 9 and P. Following Garoff and 
Meyer we obtain^ 



A/ = / - /o = 



a (T-T Q )9 2 I39 A 



99 P 2 K 
-CP9 + - EP + — \V6\ + f w , 



(4) 



where /o is the free energy density of the undistorted SmA 
phase, oq, (3 and To are the mean-field material constants, C 
is the piezoelectric coupling constant, £o is the dielectric con- 
stant, x is the electric susceptibility, K is the representative 
elastic constant of the SmC phase, and f w (6) describes the 
free energy costs at the LC-wall interface. 



Minimization of / with respect to P yields P = 
e ox(E + C9). We take this into account and introduce 

scaled dimensionless quantities s — 9^J ao ^ ^j^, Eq = 

EqC'EoX (aO(Tf-T)) ' A ^ = A / Q 2( T f-T)2' /«> = 

U a 2 {T P a _ T] 2 , V = dV, and J = f a g (r f_ r)a , where 

y _ e °x describes the second order structural phase tran- 
sition temperature for E — 0. If omitting the tildes it follows 



A/ = f c + f e + ff + f w , where 

s 2 s 4 

fe = ^ |V S | 2 , (6) 

ff = -sE cos(cJt). (7) 

The condensation term f c enforces s = So = ±1, the elastic 
term f e penalizes deviations from a spatially homogeneous 
ordering, and the term // introduces the source of the coherent 
input signal. We proceed by adopting a standard form of the 
dissipation free energy term, neglect the f w contribution, and 
after discretization we reproduce Eq. ([T), where S{ = s(fi). 
Finally, in a SSFLC cell, the dynamic Gaussian noise £j(t) 
entering Eq. (Q~|) could be produced by a randomly varying 
external electric field. 
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Next, we rewrite Eq. ([l} into a more general form 

ds 

^=8i-a$+^Ji J (s j -8i)+E+V2Dti(t)+s i Wi, (8) 

3 

where the quantities and Wi now introduce qualitatively 
different origins of static disorder. The physical origin of 
these contributions in systems of our interest is as follows. 
Standard SSFLC cells are extremely susceptible to imperfec- 
tions of the confining walls. One can substantially improve 
their mechanical stability by introducing a polymer network 
into the LC. The resulting system is referred to as the PSFLC. 
For a low enough concentration of the polymer the bistabil- 
ity is preserved. The network introduces a kind of a bulk 
field. Its local impact could be approximately modeled by the 

2 

LC-polymer coupling term f w ~ — —yS giving rise to the 
S{Wi contribution in Eq. (|8}. This term introduces, due to the 
essentially random nature of the polymer network, a certain 
degree of static randomness or disorder into the system. 68,69 
Furthermore, the presence of the polymer can give rise to 
spatially varying elastic interactions within the system, re- 
sulting in spatially dependent coupling constants Jy . Lastly, 
E = Eq cos(ui) is now a generalized weak periodic field with 
a dynamic random component approximated by the Gaussian 
noise £j(£) with the same properties as in Eq. (Q~|). 

Prior to examining the results, we further generalize Eq. ([8]l 
in order to investigate a variety of soft matter systems experi- 
encing qualitatively different origins of disorder. We consider 
influences of random dilution (RD), random bond (RB), spin 
glass (SG), and random field (RF) type of static disorder. By 
the RD type we set Jij = at randomly chosen interactions 
with a probability 7, whereas at the remaining sites it holds 
Jij = J. By the RB type we allow Jy to vary randomly 
within the interval J y - 6 [ J - A J, J], where < AJ < 2J. 
The case A J = 2 J corresponds to the SG universality class, 
where Jy exhibits random values within the interval [— J, J]. 
By the RF type of static disorder a value of the field Wi is cho- 
sen randomly within the interval [-Aw, Aw]. Furthermore, 
we consider either short-range (interactions only amongst the 
four nearest neighbors) or infinite-range (all-to-all coupling) 
interactions J y - . In the latter case we assume that a value of 
Jij is independent of the distance between the i-th and j-th 
oscillator. 



IV. RESULTS 

In what follows, we study the phenomenon of SR in the sys- 
tem described by Eq. ((S). We first focus on the influence of 
different types of static disorder via an analytical treatment of 
equilibrium values of the i-th oscillator s, within the ensem- 
ble. Afterwards, we study numerically the combined influence 
of dynamic and static disorder on the mean-field response of 
the studied soft matter system. 



A. Influence of static disorder 

We estimate how different sources of static disorder modify 
the bistable configuration enforced by the condensation term 
f c . Since bistability is a key ingredient enabling the observa- 
tion of SR, these insights will help us to understand and inter- 
pret effects of the joint impact of static and dynamic disorder 
that we are going to present in the next subsection. In order 
to assert the impact of static disorder on bistability, we treat 
disorder terms (i.e. the Jij and Wi contribution) as weak per- 
turbations. Moreover, we neglect the dynamic disorder [Gaus- 
sian noise &(£)], switch off the periodic electric field E, and 
express s, as 

Si = s + Xi, (9) 

where so = ±1. Linearization of the statical part of Eq. ([8]) 
with respect to Xi yields — 2xt + Y^j Jij( x j ~ x i) + w i( s o + 
Xi) = 0. For J^ = we obtain in the lowest approximation 

Xi ~ s Wi/2. (10) 

Accordingly, we conclude that the presence of RF type disor- 
der of any strength modifies the equilibrium configuration of 

Si. 

To estimate the influence of random variations in Jy we 
set Wi — 0. For convenience we introduce the renormalized 
random interaction coupling matrix JL as 

Jlj = Jij-5ij^2Jkj, (11) 
k 

where is the Kronecker symbol. We further express x.- L in 
the eigenbasi a 70 ' 71 of the matrix Jy 

Xi=J2 A ^\ (12) 

a 

(13) 

3 

where £! Q are the eigenvalues, ef are the corresponding 
eigenvectors and A a are the weight constants. The eigenvec- 
tors are normalized according to 

i 

We expand the statical part of Eq. ([8]) for E = including 
the cubic term in Xi, take into account Eqs. ( fT2b and ( fl~3b . 
multiply the resulting equation with e\ , and sum it over the 
lattice sites. It follows 
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(15) 



Q,/3 



r 



where 



s («)JfflJ>l) 



(16) 



Neglecting the coupling among different eigenmodes and 

(ri) 

assuming — one obtains the equation for the 77-th mode 
amplitude 

A*F[ v) + Ar,(Q - 17„) = 0, (17) 
where 17 = 2. It follows A n (£l v < 17) = and 



^,,(17,, > 17) 



i 7 



('/) 



(18) 



The corresponding effective Hamiltonian H v determining the 
?y-th eigenmode amplitude can be defined as 



fir. 



4 4 



(19) 



where 17 = 2 plays the role of the generalized temperature. 
On decreasing 17 the eigenmodes open via a continuous phase 
transition at 17 = 17^ and below 17^ its amplitude is given by 
Eq. ( TT8T >. It follows that the disorder in can trigger struc- 
tural changes only for large enough disorder strengths, i.e. at 
least one eigenvalue must be larger than 17. When departing 
from the bistable solution preferred by f c only modes charac- 
terized by 17^ > 17 are thus present. 

Note that in the derivation of Eq. ( fT8l we have assumed that 
the eigenmodes do not overlap. This is partly justified only for 
localized modes 71 . These are expected for short-range inter- 
actions in Jij and 17, ; > 0. On average they appear at different 
regions and consequently avoid overlapping. The response of 
the system is in this case given by 



Si = s 



A f ( 



(a) 



(20) 



where only the modes with 17 a > 17 = 2 participate. The 
extended modes thus appear only for 17 r) ~ 0. 72 

In case of long-range interactions in Jy the modes are 
extendedjS, For example, for infinite-range interactions and 
SG type disorder the distribution of eigenvalues stays within 
the interval 17 ?) € [— 17 max , 17 max ] according to the symmetric 
Wigner probability distribution^ 2 . 



7rl7 n 



17 2 - 17 2 

J 'max J '77 > 



then the modes would appear via a discontinuous transitions 
as a function of 17. 

To illustrate the influence of static disorder and different 
interaction ranges on Xi configurations we focus on the quan- 
tity l/F^\ where F^ v) is defined by Eq. <Q2]). It estimates 
the number N v of sites participating in the 77-th mode, i.e. 



(</) 



N v . The estimate follows from the normalization 



condition given by Eq. ( fT4b . We assume that all 'open' sites, 
characterized by x^ 7^ 0, are comparable in magnitude. 
Therefore {x^l ~ 1/a if site i is opened, whereby a stands 

for a positive constant, but otherwise x^ = 0. In Eq. ( TBI the 
sum is restricted over N v open sites, resulting in N v /a 2 = 1. 



'N n and F 4 



(n) 



E ( (e<*>) 



Consequently, 

1/N V . Non-localized (extended) modes are characterized by 

N„ 



N and localized by N v < N. 



In Fig.|2]we plot l/F^ ~ N v as a function of mode eigen- 
values Cl v for N — 1000. We consider the RB case with 

A J = J, and the SG limit A J = 2 J. The resulting l/F^ v) 
values are shown for short- and infinite-range interactions in 
the top and bottom panels, respectively. By short-range in- 
teractions the modes are typically localized [in Figs. [2 a) and 
(b) N v is typically less than 10]. On the contrary, by infinite- 
range interactions most states are extended [characterized by 
N v > 100, as shown in Figs. [2fc) and (d)]. We further find 
that by increasing the value of A J, where Jy e [ J — A J, J], 
the maximal value 17 max increases. According to the estimate 
given by Eq. ( fT6b the modes can open if 17 max > 17 = 2. 
Our calculations suggest that such modes are enabled only for 
A J > J. Consequently, we find that the disorder in no- 
tably affects the SR phenomenon only if A J > J. Note that 
Eq. (IT~6b holds only for relatively weak static disorders via 
Jij, and for the case that modes do not strongly overlap. 

We further show that the disorder introduced via the cou- 
pling interaction Jy is unlikely to destroy the bistability en- 
forced by the condensation term. For this purpose we origi- 
nate from Eq. (0 and set E — 0, Wi — 0, and = 0. We 
express Sj as 



S + 6si 



(22) 



where Ssi describes fluctuations around the mean field S of 
the system. If the bistability is preserved, then S = ±|5|. 
In the opposite case S = 0. We sum the resulting equa- 
tion over all lattice sites, assume -i < ^ff 1 ~ j? J2i ^Si ~ 



0, and obtain the equation 



where 17 max stands for the maximal eigenvalue. Therefore, a 
single mode is expected to dominate in Xi, and this mode hin- 
ders the opening of the remaining modes. Further, if F^ 7^ 



S - S 3 - 3SM = 0, (23) 
where M measures the strength of departures about the aver- 



6 




age: 

M =^E^- w 

i 

It follows that the bistability is preserved for M < 1/ 3, where 
S = ±Vl - 3M. (25) 

Finally, within this subsection we address the condition 
where the bistability is lost with eigenmodes of the matrix 
JL. From Eqs. d22i > and (O it follows Ssi = sq — S + Xi. 
The bistability is lost when M = 1/3 and S = 0. With this 
in mind, using the definition given by Eq. (l24l and expansion 
given by Eq. (TT2l . we obtain the condition 

mm m 

(26) 

where the sum runs over the open modes, the amplitudes of 
which are estimated by Eq. (TT~8T >. It is obvious that the condi- 
tion given by Eq. d26| i cannot be satisfied. Therefore, the static 
disorder in the coupling strength cannot, on its own, destroy 
the bistability. 



B. Combined influence of static and dynamic disorder 

In this subsection we analyze numerically the response of 
the system in terms of the mean field S if both static and dy- 
namic disorder are introduced. We use periodic boundary con- 
ditions for our simulations and consider short- and infinite- 
range interactions in Jy . In the former case only the nearest 
neighbors are coupled, whereas in the latter the coupling is 
independent of the physical distance between the units. By 
considering both short- and infinite-range couplings we take 
into account two limiting cases of interactions that interpolate 
between the 3D structures usually entailed within soft matter 
systems. Importantly though, we found the actual usage of 3D 
models prohibitive due to substantial computer resources that 
would be needed to simulate them. We do, however, demon- 
strate that our results are independent of the system size in 
Fig. [6] and thus should be readily observed also in the ther- 
modynamic limit. 

We start by examining the temporal evolution of the mean- 
field response S for three different levels of Gaussian noise 
D. As the reference state we consider a configuration with 
Jij = J = 4 in the absence of any additional static disor- 
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FIG. 3: Temporal evolution of the mean-field response S of the sys- 
tem by different levels of (dynamic) Gaussian noise D. (a) D — 0.7, 
(b) D — 0.9 and (c) D = 1.3. Gray lines depict the periodic driving 
(amplitude is scaled to 1). In panel (b) the resonance condition is ful- 
filled, resulting in the optimal correlation between the weak periodic 
driving and the response of the system. Employed system parameters 
are: JV, = J = 4, w t = 0, E = 0.04, w = 0.01 and N = 100. 



der and infinite-range of interactions. Figure [3] features the 
results. In panel (a) the noise level is clearly to weak to as- 
sist the periodic field strong enough to induce flips between 
the two minima of the potential (note that Eq = 0.04 assures 
that E itself is too weak to induce the flips). In sharp contrast, 
panel (c) depicts a heavily noisy response, lacking a clear tem- 
poral structure that would indicate any particular correlations 
with the weak deterministic driving. Only in panel (b) the cor- 
relation between S and the weak periodic field E is evident 
already with the naked eye, and indeed, visually manifests the 
SR phenomenon in the examined soft matter system. To sup- 
port this visual assessment, we examine S evoked by different 
D via the Fourier coefficients Q [see Eq. ([3])] that quantify the 
correlation between S and E. Figure [4] depicts a typical bell- 
shaped dependence of Q, thus confirming that an intermediate 
level of Gaussian noise is able to optimally assist the weak pe- 
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FIG. 4: SR in the examined system. By an intermediate value of D 
(denoted by the dashed vertical line) the correlation of the mean-field 
response S with the weak periodic driving E is maximal, resulting 
in a bell-shaped dependence of Q. System parameters were the same 
as in Fig. [5] 



riodic driving E to induce flips between the two minima of the 
potential. This is the hallmark of the SR effect, which we thus 
demonstrate in the examined soft matter system. 

In order to distinguish impacts of different types of static 
disorder on the SR resonance effect exemplified in Fig.|4j we 
continue with examining the mean field responses S via the 
Fourier coefficients Q. First, we consider the influence of de- 
creasing coupling strengths Jy = J to obtain a reference for 
further simulations. In particular, note that both the RD and 
RB types of static disorder progressively lower the average 
coupling strength as their magnitude increases, and thus 
the reference simulation in which the interaction strength be- 
tween all the coupled oscillators J is gradually reduced may 
provide important pointers for the understanding of the effect 
of RD and RB static disorder. The resulting responses quanti- 
fied via Q in dependence on D and J are shown in Fig. |5ja) 
separately for short- and infinite-range interactions. For,/ = 4 
the correlation of S with E peaks at D = D r sw 0.35 (short- 
range) and D r « 0.9 (infinite-range; see also Fig. |4j, whereby 
the subscript r marks the resonance value of D. In both cases 
the resonance peak shifts toward lower values of D — D r as 
the coupling strength is reduced, and moreover, an accompa- 
nying increase in peak values of Q prior to reaching J = can 
be inferred in case of infinite-range interactions (the global 
maximum of Q is obtained at D = D r = 0.11 and J = 0.17). 
At a glance, it can be observed that RD [Fig. HJb)] and RB 
[Fig. |3c)] types of static disorder by infinite-range interac- 
tions (right column) deliver qualitatively identical results as 
obtained by the reference simulation [Fig. [2 a)]. Indeed, by 
increasing 7 (RD) or A J (RB) the average coupling strength 
Jij decreases, which due to the fact that by infinite-range in- 
teractions each unit experiences approximately the same cou- 
pling with others, simply results in the same effect as if J 
would be gradually reduced equally for all units. Accordingly, 
in the right (b) and (c) panels the optimal D — D r decreases 
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FIG. 5: Response of the system, quantified via Q [see Eq. y]l], in dependence on the Gaussian noise level D as well as the type and strength 
of static disorder. Results are presented separately for short-range (left column) and infinite-range interactions (right column). Rows depict: 
(a) impact of different Jij = J; (b) RD; (c) RB; and (d) RF. In all panels the color profile is RGB linear (blue > green > yellow > red), blue 
depicting 0.0 and red 0.45 values of Q. The employed system size was TV = 100 and Q were evaluated over n — 1000 oscillation periods of 
E. In addition, all results were averaged over 50 different realization of every type of static disorder. Where applicable, other parameters were 
the same as in Fig. [3] 



as 7 and A J increase, and global maxima (across the whole 
2D parameter space) of Q are obtained prior to reaching the 
limits of the RD and RB universality class. Specifically, the 
global maxima of Q are obtained at D = D r = 0.11 and 
7 = 0.96 by the RD case, and by D = D r = 0.18 and 
A J = 7.1 by the RB case. Conversely, by short-range inter- 
actions the individuality of couplings amongst units is more 
pronounced, and accordingly, the reference simulation in the 
left column of Fig. [5ja) bears less similarity with the impact 
triggered by RD and RB types of static disorder, especially if 
compared to infinite-range interactions. In case of short-range 
interactions the RD and RB types of static disorder introduce 
effects that go beyond a simple decrease in J. The main trend, 
however, remains the same in that the resonance peak of Q 
shifts towards lower values of D — D r as the magnitude of 
static disorder increases, but nevertheless, the global maxima 
of Q are obtained by 7 = and A J = 0. Moreover, sub- 



stantial removal of links with probability 7 by the RD case 
seems to have a rather profoundly deteriorating effect on the 
maximally attainable peak value of Q. 

Most notable, however, is the qualitatively different re- 
sponse that is evoked by the RF type of static disorder de- 
picted in Fig. EJd), which sets in irrespective of the range of 
interactions since Wi always has only a localized impact on 
each particular unit. It can be inferred that by the RF type 
of static disorder the resonance peak shifts toward higher val- 
ues of D — D r as Aw increases, and also, an accompanying 
decrease in peak values of Q can be inferred (global max- 
ima of Q are obtained at Aw = 0). This is significantly 
different from what we have observed by RD and RB types 
of static disorder, especially so by infinite-range interactions. 
The discrepancy can be explained by our preceding separate 
treatment of effects of static disorder. In particular, unlike 
variations in disorder introduced via Jy does not change 
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AJ=8 



M=J 




for the full range of A J are presented in Fig. |6ja). Shaded 
are areas in the D — A J parameter space where Q > 0.18. 
Clearly, the SR is preserved as N increases. To enable a bet- 
ter quantitative insight of the preservation of SR, we show in 
Fig. |6jb) excerpts of the resonant behavior for A J = J. In- 
deed, it can be observed nicely that not only the combined 
effects of static and dynamic disorder are preserved, but also 
that the resonance behavior saturates for large enough N. This 
is in accordance with previously reported results of Monte 
Carlo simulations on a simple Ising-type lattice model^i thus 
fully confirming the generality of our findings with respect to 
variations in system size. 



V. SUMMARY 




0.0 0.2 0.4 0.6 o.; 

D 



1.0 1.2 1.4 



FIG. 6: loint effects of RB type static disorder and Gaussian noise in 
dependence on the system size N by infinite-range interactions, (a) 
Shaded outlines, representing Q > 0.18, of the resonance behavior 
for different N. (b) Excerpts of SR depicted in panel (a) obtained by 
A J = J. Importantly, the SR is preserved and saturates in magni- 
tude above N ~ 300. 



the double well potential of Sj if s$ = Sj since such pairs are 
preferred for a positive value of Jy , and hence the essential 
bistability property for SR remains intact. On the other hand, 
disorder introduced via Wi directly influences the energy land- 
scape experienced by Sj, and thus destroys bistability for large 
enough Aw. This in turn explains the ever increasing values 
of D = D r needed for a resonant response and the deteriorat- 
ing peak values of Q depicted in both panels of Fig. [2d) as 
Aw increases. 

Finally, we examine whether above results are qualitatively 
independent of the system size. The latter is an important 
issue since by soft matter systems the thermodynamic limit 
N — > oo can in general be considered fulfilled. For the pur- 
pose of this analysis we consider the RB type of static disorder 
as a representative example (for other types of static disorder 
identical results are obtained), and moreover, focus only on 
infinite-range interactions as variations in system size do not 
influence the SR in case of short-range interactions. Results 



We studied the phenomenon of SR via variations of static 
and dynamic disorder in a soft matter system. As a representa- 
tive example, we choose a polymer stabilized ferroelectric LC 
cell, which is an essential ingredient in several electro-optic 
applications. The cell's thickness was considered shorter than 
the pitch of the helix, which resulted in a bistable dynamics 
governing the discretized equations of the LC configuration in 
the SmC phase with qualitatively identical properties as exhib- 
ited by the paradigmatic bistable overdamped oscillator. As 
the latter is a generic model for SR, the observation of noise- 
enhanced correlation between a weak periodic field and the 
PSFLC response could thus be expected. Indeed, by apply- 
ing a subthreshold electric field E along the surface normal of 
the two plates that alternatively favored the competing Up and 
Down configurations, we were able to demonstrate the phe- 
nomenon of SR in a typical soft matter system by measuring 
its mean-field response S. 

More precisely, we focused on the combined influence of 
static and dynamic disorder, whereby the former arises natu- 
rally due to the random character of the polymer network that 
is introduced into the LC in order to improve its mechanical 
stability. In particular, the coupling between LC molecules 
can then no longer be considered uniform because of induced 
local variations in the effective temperature within the system, 
thus introducing RD and RB types of static disorder, and be- 
cause of the branched polymer structure introducing RF type 
disorder. Importantly, we generalized the properties of the ex- 
amined LC in order to obtain broader insights into the behav- 
ior of qualitatively different systems, which are however, gov- 
erned by the same generic equations. In addition to consider- 
ing these three types of static disorder, we also examined the 
system's dynamics in dependence on the system size as well 
as on local- and infinite-range interactions. First, we showed 
analytically that static disorders introduced via the coupling 
strength (RD and RB) and via the random field Wi (RF) 
affect the system's bistability in qualitatively different ways. 
While the RF static disorder is very effective in destroying the 
bistable character of the system even by small magnitudes, 
the RD and RB types of static disorder fail to have effects un- 
less the level of randomness is large. We further show that 
the range of interactions strongly affects the response of the 
system with RD or RB type of disorder. These analytical es- 
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timates were obtained by taking into account a system with 
no dynamic disorder and relatively weak static disorder. Our 
numerical simulations, where the combined influence of static 
and dynamic disorder was taken into account, confirmed our 
analytical estimates. While RD and RB types of static disor- 
der can be simplified mostly to decreasing coupling strengths 
Jij = J (with some differences emerging with respect to the 
range of interactions amongst LC molecules), the RF type dis- 
order induces a qualitatively different response. Particularly, 
RD and RB (becomes SG in the maximal disorder limit) de- 
crease the level of dynamic disorder warranting the optimal re- 
sponse, whereas the RF evokes exactly the opposite effect by 
increasing the optimal D that is needed to resonantly fine-tune 
S in accordance with E. These observations are shown to be 
independent of the system size and largely independent also 



of the range of interactions, which implies that they are gen- 
erally valid and potentially applicable beyond the presently 
employed setup. Indeed, our study leads us to conclude that 
soft matter systems are viable candidates for additional theo- 
retical as well as experimental research, potentially leading to 
sensitive detectors with fortified capabilities due to the phe- 
nomenon of SR. 
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